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Abstract 

We review the description of nucleon structure functions in the instanton vacuum. 
This includes the calculation of the twist— 2 parton distributions at a low normaliza- 
tion point as well as higher-twist matrix elements. The instanton vacuum with its in- 
herent small parameter, the packing fraction of the instanton medium, p/R, provides 
a consistent picture of the non-perturbative gluon degrees of freedom at the scale 
~ 600 MeV. The twist-2 quark and antiquark distribution are of order unity, 
while the twist-2 gluon distribution is of order (p/i?) 4 . Twist-4 matrix elements de- 
termining power corrections to the Bjorken, Ellis-Jaffe and Gross-Llewellyn-Smith 
sum rules are found to be of order (p/R)°. We present numerical estimates for the 
parametrically large quantities. 



f Talks presented at the XXXVII Cracow School of Theoretical Physics, Zakopane, May 30 - June 10, 
1997. 



1 



In these talks we give a brief summary of recent progress in understanding the deep- 
inelastic structure of the nucleon, both at leading and non-leading twist level, in the 
instanton vacuum [I], |2|, |3], |J. Our aim is to show that the instanton vacuum, with its 
inherent small parameter — the packing fraction of the instanton medium, p/R — provides 
a basis for a consistent and quantitative description of nucleon structure functions. 

Leading and non-leading twist. The non-perturbative information which enters in the 
QCD description of deep-inelastic scattering and other related experiments is contained 
in nucleon matrix elements of operators of twist 2 and higher. The moments of the non- 
power suppressed part of the structure functions are given by matrix elements of operators 
of leading twist; in the unpolarized case these are the twist-2 quark and gluon operators 
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and similarly for the polarized case, see ||. Here, t$ = 1,tns — t 3 are flavor matrices. 
Alternatively, one may work with non-local (light-cone) operators, which serve as gener- 
ating functions of the series of local twist-2 operators. In a partonic language, the matrix 
elements of these operators can directly be interpreted as parton distribution functions^, 
the scale dependence of which is described by the DGLAP evolution equation: 
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where n denotes a light-like four-vector, n 2 = 0. Operators of higher twist arise in the 
description of power corrections f|. For example, in the unpolarized case the l/Q 2 -power 
corrections to the Gross-Llewellyn-Smith sum rule are governed by the matrix element of 



^^Here, q_f(x) corresponds to the quark distribution at positive x, and to minus the antiquark distribution 
at negative x. 
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the twist -4, spin-1 operator (Mjv is the nucleon mass) 



lJ2(P\^ al5 F^\P) = 2M 2 N C { i ) P^. (5) 



spin 



In polarized scattering the l/Q 2 -power corrections to the isovector and isosinglet com- 
binations of the first moment of the polarized structure function g\ — the Bjorken and 
Ellis- Jaffe sum rules — involve the matrix elements of the twist -3, spin-2 operators || 

(PS\$ t ns ,s {l a F Pl + ^F^) i[>\PS) - traces (6) 
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(the same matrix elements contribute also at leading twist level to the third moment of 
the structure function g 2 ), and the matrix elements of the twist-4, spin-1 operators, 

(PS\^r NS , sla F^\PS) = 2M 3 N f£{ s SP. (7) 

Here, S is the nucleon polarization vector, S 2 = —1. The operators here are assumed to 
be normalized at scale //; the scale dependence of the matrix elements is described by the 
renormalization group equation. The twist— 2 quark-, antiquark- and gluon distributions 
at a low normalization point have been determined by fits to data from a variety of 
experiments The twist— 3 matrix element has recently been extracted from 

measurements of the structure function g% ||. Experimental knowledge of the higher-twist 
matrix elements entering only in power corrections, such as the twist-4 matrix element 
f( 2 \ is still rather poor |fL0| . 

Any attempt to calculate the matrix elements mentioned here from first principles re- 
quires an understanding of the non-perturbative effects giving rise to the structure of the 
nucleon. While the gross features of the leading-twist quark distributions can be under- 
stood in phenomenological models like the quark model or the bag model, to describe the 
gluon distribution or higher-twist matrix elements explicitly involving the gluon field one 
needs a theory of the non-perturbative fluctuations of the gluon field. Lattice calculations 
of structure functions have been making steady progress during the last years; however, 
they are still far from giving a satisfactory quantitative description JTT] . 



Instanton vacuum. A microscopic picture of the non-perturbative fluctuations of the 
gluon field is provided by the instanton vacuum. For an introduction to the instanton 
vacuum and its applications to hadronic physics we refer to the extensive literature on this 



subject, e.g. the recent reviews [[12], |13j; we can here touch upon only those aspects directly 
relevant to structure functions. Instantons and antiinstantons (/ and / for short) are 
particular field configurations which are solutions of the Euclidean Yang-Mills equations, 
characterized by a size, p, center, z and color orientation given by an SU(N C ) matrix, U, 

Al{x; z,U)i { l) = W-z)O ab (^)l, 

Ux) = {x2 2 f p2)x2 ^ 0« b = Itr[A^A b Wt]. (8) 
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Here, {jj T ) b ^ u = fj b , are the 't Hooft symbols. Instantons have many special properties; 
not all of them are of interest to us here. Let us note only that a single 1(1) is an 
0(4)-symmetric field configuration — this fact makes for important differences between 
instanton contributions to operators of different spin, see below. 

In the instanton vacuum one considers non-perturbative effects due to field configura- 
tions with a finite density of 7's and J's. The medium of J's and J's stabilizes itself due 
to instanton interactions |]15|| , meaning that the average size of the instantons in medium 
is finite, 

p ~ (600 MeV)- 1 . (9) 

The coupling constant is fixed at a scale of order p~ l , so when we evaluate matrix elements 
of QCD operators below it is implied that the operators are normalized at p ~ p _1 . It 
should be stressed that no external scale is introduced here; all parameters of the instanton 
medium are obtained in terms of the QCD scale parameter, Aqcd- Hence this approach 
preserves the essential renormalization properties of QCD. 

The most important property of the instanton vacuum is the small packing fraction of 
the medium, i.e., the small ratio of the average size of the instantons in the medium to the 
average separation between nearest neighbors, p/R ~ 1/3 |TJ|, |T5] . This small parameter is 
the starting point for a systematic analysis of non-perturbative phenomena in this picture. 

In particular, the instanton vacuum explains the dynamical breaking of chiral symme- 
try. The Dirac operator in the background of one 1(1) has a localized zero mode, 

" 3 + fi I{I) (x; z,U)] $±(x; z,U) = 0. (10) 



In the medium the zero modes associated with the individual instantons delocalize |16| 



resulting in a finite fermion spectral density at zero eigenvalue, which by the Banks- 
Casher theorem is equivalent to the chiral condensate [FT]]. Alternatively, one can derive 
the effective action of fermions in the instanton medium in the l/A^-expansion, integrating 
over the instanton coordinates in the ensemble [IB]. An individual 1(1) interacts with the 



fermion field through the zero mode, i.e., through a "potential" 

V I{ l0f,i>f] = Jd 4 xJd 4 y^ f (x)^ ± (x;z,U)^ ± (y;z,U)^ f (y). 
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In leading order in p/R the effective action exhibits a many-fermionic interaction, which is 
given by the one-instanton average of eq.(|Tl~P and has the form of the 't Hooft determinant 
in flavor indices [|H| , 

N f i _^ 

J d A zdU\{V I{I) [i, f ,^ f } cc det^—J^. (12) 

In addition, there is a form factor (not written here for brevity) related to the finite size of 
the zero-mode wave function, which makes the interaction vanish for quark momenta larger 
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than p _1 . Chiral symmetry is spontaneously broken due to this many-fermionic interac- 
tion: the quarks acquire a dynamical (momentum-dependent) mass, a quark condensate 
develops, and a massless pion appears as a collective excitation. In the l/iV c -expansion 
one can easily construct the effective action in the chirally broken phase. It can be formu- 
lated as a theory of massive quarks, interacting with the pion field in a chirally invariant 
way, 

Z = jVir J Vi)V^) exp / d 4 x i^(x) [i$ + iMe^^] ij;(x). 

(13) 

The effective theory applies for quark momenta up to the inverse instanton size, pr l ~ 
600 MeV, which acts as a cutoff. For the applicability of this effective theory it is crucial 
that the ratio of the dynamically generated quark mass to the cutoff is proportional to the 
packing fraction of the instanton medium, 

Mp oc (|) 2 . (14) 

Hence the diluteness of the instanton medium guarantees that the picture of massive 
"constituent" quarks applies in a parametrically wide range of momenta. Finally, we note 
that the nucleon is obtained in the l/iV c -expansion as a chiral soliton of the effective 
theory, eq.(Q) 50] . This picture of the nucleon gives a very reasonable description of 



a variety of hadronic properties such as the N—A splitting, electromagnetic formfactors, 
axial coupling constants etc. [pifl . 

The instanton vacuum, with the resulting effective chiral theory, allows to evaluate 
hadronic matrix elements of QCD operators involving the gluon field. In ref. |1|] a method 
was developed by which "gluonic" operators can systematically be represented as effective 
operators in the effective chiral theory. Let J 7 [A] be a gluonic operator, i.e., some function 
of the gauge field. To leading order in the packing fraction, p/R, the interaction of the 
gluon operator with the fermion field is mediated by single instantons. The effective 
operator, denoted by "JF", is obtained by substituting in J-[A] the gauge field of one 1(1) 
and integrating over the collective coordinates, 

Nf 

".F'[<M oc ]T [<PzdUF[A IiT) (z,U)] l[V I{I f ^ f }. (15) 

I+I f 

In higher orders of p/R one needs to take into account many-instanton contributions to 
the effective operator. Below we shall need the effective operator corresponding to the 
gauge field itself, which is given by (cf. eq.(]^)) [|J 

"^"W oc J2 jd*zf u {x-z)4>{z)^a^-^e ± ^ a ^{z). 

(16) 
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Again we have suppressed the form factors coming from the zero mode wave function. 
Note the presence of the pion field, as a result of which the effective operator eq.(|T6"l) is 
chirally invariant — as it should be, since the gluon field is flavor neutral. 

It is important to note that the representation of QCD operators as effective operators 
is possible relying entirely on the approximations already inherent in the effective theory - 
the diluteness of the instanton medium and the l/iV c -expansion; no additional assumptions 
are required. It was shown in |fj] that this approach preserves the essential renormalization 
properties of QCD; for example, the QCD trace and U(l) anomalies are realized at the 
level of hadronic matrix elements. This method is thus well suited for computing matrix 
elements of the QCD operators of twist 2 and higher twist of interest here. 

Twist-2 matrix elements in the instanton vacuum. Let us first consider the twist-2 
gluon operators, eq.(^). Since we are interested only in forward matrix elements we may 
average the operator position over the 4-volume. The definition of the effective operator, 
eq. (|i5|) , implies the integral over the instanton coordinate, z. Since the instanton field is 
0(4) symmetric the only tensor one has at hands to construct the effective operator is 
the Kronecker delta, but from it is impossible to construct a traceless symmetric tensor! 
(When working with light-like components this follows from 5 ++ = 0.) Thus, the effective 
operators for the twist-2 gluon operators, eq.(@), or, more generally, the nonlocal operator, 
eq.(H), vanish at one-instanton level. We conclude that the twist-2 gluon distribution is 
parametrically suppressed in p/R. 

To determine the twist-2 gluon distribution quantitatively one has to include (at least) 
the two-instanton contribution to the effective operator. Preliminary results indicate that 
the gluon distribution is, in fact, proportional to (Mp) 2 , which is parametrically of order 
(p/i?) 4 .0 It is interesting to note that, numerically, a suppression of the gluon relative to 
the singlet quark momentum fraction by a factor (Mp) 2 ~ 0.3 is consistent with the GRV 
parametrization of the data at a normalization point of p ~ 600 MeV ||. 

In the twist-2 quark operators, eq.(P, the gauge field enters through the covariant 
derivative, or, equivalently, through the path-ordered exponential in the non-local opera- 
tor, eq.@. It is interesting to ask how much the gauge field contributes to the moments 
of the quark distribution functions in the instanton vacuum, which is formulated in the 
so-called singular gauge in which the instanton field has the form eq.@. For simplicity, let 

(2) 

us consider the second moment of the singlet unpolarized quark distribution, A s , which 
is given by the matrix element of the operator 

i(P\ipj {lil (d^ 2} -i^Al 2} U\P) -traces = 2^ P m - traces. 

(17) 

It is instructive to compute the matrix element not immediately in the nucleon, but first in 
a "constituent" quark, i.e., the massive quark of the effective chiral theory. It is easy to see 

2 The gluon distribution in the instanton vacuum has recently been studied in a different approach by 
Kochelev p3]. However, the contributions taken into account there do not represent the full answer to 
order (p/Kp. 
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that the short derivative in eq. (|T 



(2) 

makes an order unity contribution to A s . Computing 
the gauge field contribution to the matrix element using the effective operator, eq. flloD , 
one finds that it is of order (Mp) 2 
the short derivative. Thus we have 



oc (p/R) , i.e., parametrically suppressed relative to 
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This is consistent with the fact that the gluon distribution is of higher order in the packing 
fraction: To order (p/R)° the quarks carry the entire momentum, and the gluon distribu- 
tion is zero. (To compute the full {p/R) 4 contribution would, again, require to take into 
account two-instanton contributions to the effective operators, as well as to the effective 
quark action.) The above statements are easily generalized to higher moments. In fact, 
speaking of parton distribution functions "inside the constituent quark" one may say that 



^2 <7/( X ) quark 
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9v^/ quark 



8(x - 1) + O 
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To leading order in p/R the constituent quark has no structure. Consequently, in lead- 
ing order in p/R it is justified to identify the distribution of "constituent" quarks and 
antiquarks in the nucleon with the actual parton distribution at the scale p ~ p . In 
another way of saying, when computing the twist-2 quark distribution function in the 
effective theory in leading order in p/R one can identify the QCD quark fields normalized 
at p ~ pr l with the quark fields of the effective chiral theory and put the path-ordered 
exponential in eq.(|3]) to unity. This is the "quarks-antiquarks only" approximation which 
was employed to compute the quark and antiquark distributions of the nucleon in the chi- 
ral quark soliton model 0. We have thus seen that this approximation has a parametric 
justification in the instanton vacuum. 

The twist-2 quark and antiquark distributions of the nucleon computed in the chiral 
quark soliton model satisfy all general requirements, such as positivity, proper normaliza- 
tion etc. This is a fully field-theoretic description of the nucleon, which, in particular, 
makes possible a consistent calculation of the antiquark distributions. All partonic sum 
rules (baryon number, isospin, momentum, Bjorken sum rule) are satisfied within the 
model. Computed so far were the leading distribution functions in the large- N c limit, the 
isosinglet unpolarized and isovector polarized ||, as well as the isovector transverse po- 
larized distribution 0]. The isosinglet unpolarized quark and antiquark distributions are 
shown in Fig.|l|, together with the GRV parametrization || ; for the polarized distributions 
see the original papers. 

To order (p/R)° the constituent quarks have no structure, and the constituent quarks 
and antiquarks carry the entire nucleon momentum. As a result the singlet quark distri- 
bution calculated in this approximation is generally larger than the GRV parametrization, 
which includes gluons at the low normalization point, see Fig.l. In higher orders of p/R 
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Figure 1: Solid lines: The isosinglet unpolarized valence quark and antiquark distri- 
butions computed in the chiral soliton model of the nucleon ||. Dashed lines: NLO- 
parametrization of GRV (p 2 = 0.34 GeV 2 ) @. 

one starts to systematically resolve the structure of the constituent quark in terms of the 
original QCD degrees of freedom, and the nucleon momentum gets distributed among 
quarks and gluons at the low scale. A 30% gluon momentum fraction at the low scale 
obtained by GRV || is consistent with the gluon distribution in the instanton vacuum 
being suppressed by a factor (Mp) 2 ~ 0.3. 

To summarize the discussion of twist-2 operators, one may say that at twist-2 level 
the effects of the instanton medium are essentially contained in the dynamical quark 
mass generated in the dynamical breaking of chiral symmetry. Instanton contributions 
to the twist-2 operators, in the sense of effective operators, eq. (|T5|) , are parametrically 
suppressed. This is what one could call a "constituent" quark picture. We note that, 
contrary to other approaches where the "constituent" quark is a largely philosophical 
object, in the instanton vacuum this term has a well-defined meaning, thanks to the 
parameter p/R P3"[ . 

Higher-twist matrix elements in the instanton vacuum. From the above discussion one 
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may have the impression that in the instanton vacuum "gluonic" contributions to opera- 
tors are always parametrically suppressed relative to quark operators. This is not so - 
instantons can make order (p/R)° in operators of higher twist. The most immediate way 
to convince oneself of this is to consider a particular higher-twist operator whose matrix 
elements vanishes by the QCD equations of motion. In addition, this exercise provides a 
beautiful check for the consistency of the effective operator method. Consider the twist-4 
matrix element obtained by projecting the operator in eq.([T7|) not on spin two but on spin 
zero: 

i(P\^Y^\P) = (P\i>{i$ + fi)ip\P) = 0. (21) 
In QCD the matrix element is zero by virtue of the QCD equations of motion. In the 



instanton vacuum we find, computing the contribution of the gauge field in eq.(|21|) using 
the effective operator, eq. fllBf ), one may show that 

(P\$(i$ + "fi")i)\P) = (P\4>(i$ +iMe il5T ^ a{x) )i)\P) = 0. (22) 

(The calculation is actually rather involved; see |§] for details). Here the gauge field 
contribution to the operator is of order unity, i.e., of the same order as that of the short 
derivative. As a result the QCD operator of eq.(pl|) reduces to an effective operator which 
vanishes identically due to the equations of motion of the effective chiral theory. From 
this we learn two things: First, instantons can make order unity contributions in twist-4 
operators. Second, the method of effective operators preserves a principal feature of QCD: 
matrix elements of operators which are zero in QCD due to the QCD equations of motion 
are automatically zero in the effective theory. We note that also other operators whose 
forward matrix elements vanish in QCD, 

(P\i; la F^\P), (PS\^ al5 F^\PS), (23) 

give zero matrix elements when translated to the effective chiral theory. 

We can now turn to the calculation of the matrix elements of twist -3 and 4 operators 
appearing in power corrections, eqs.(|, || [7]). Again, the qualitative features can be seen 
by studying the matrix elements in "constituent" quark states. Computing the matrix 
elements of the effective operators corresponding to eqs.(|5|, |7|) one finds that twist-3 
and 4 matrix elements are of different order in the packing fraction: 

twist 4: fW,C^ ~ (Mpf ~ (|Y, 

twist 3: ~ (MpflogMp ~ (f^ log (f) • 

Again, the reason for this can be seen in the 0(4) symmetry of the individual /(/). Note 
that the parametric order of the matrix elements is determined by two factors: i) the 
number of instantons participating in the effective operator (here we have included only 
the one-instanton contribution), and ii) the dependence of the quark loop integrals in 
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Instanton vacuum [|J 
Sum rules 



24] 



Sum rules p5fl 
Bag model 12 
Lattice [11] 
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Sum rules 
E142, E143, E154 
Power corr. HlO 



J NS 

-0.10 
-0.20 
-0.072 
0.11 



0.10 ±0.28 



d 



w 

NS 



-3 



~ 10 

0.072 
0.072 
0.063 
-0.13 

0.043 ± 0.046 



C 



w 



0.36 



0.37 



scale/ GeV 2 
~ 0.4 
1 
1 
5 
4 
1 
3 
1 



Table 1: Numerical results for the flavor-nonsinglet spin-dependent twist-4 and 3 matrix 
elements fj^g and d NS , eqs.(|| |7|), and the flavor-singlet spin-independent twist-4 matrix 
element Cg , eq.(|^). Shown are the results obtained from the instanton vacuum ||, 



from QCD sum rule calculations [24| |27|, from the bag model p6 |, and from lattice 



calculations ITTH . Also shown are estimates of d$ s based on measurements of the structure 

p(2) 



function g 2 ||, and estimates of f^' s from an analysis of power corrections to g\ |T0 



the matrix element (obtained by closing the quark lines on the many-fermionic effective 
operator) on the cutoff, p" 1 , keeping in mind that (Mp) 2 ~ (p/i?) 4 . An order unity 
contribution can come only from integrals which are "quadratically divergent", meaning 
they are proportional to p~ 2 . Incidentally, this last fact implies that, in the framework 
of our effective theory, the dominant contributions to higher-twist matrix elements come 
from "divergent" loop diagrams where the effective many-fermionic operator couples to 
a single constituent quark, not from diagrams describing interactions of more than one 
constituent quark mediated by the many-fermionic effective operator. In the latter all 
momenta are cut by the bound-state wave function of the nucleon, not by the cutoff, 
fT l . Losely speaking, one may thus say that in our picture higher-twist matrix elements 
measure properties of the individual constituent quarks, not correlations between them. 

We have computed the leading higher-twist nucleon matrix elements in the large--/V c 
limit; see M for details. These are the flavor nonsinglet spin-dependent ones (d^ s , fffi s ) 
and flavor singlet spin-independent ones (C ( s>). Results are shown in table IB. The result 

(2) 

for d NS should be taken as an order-of-magnitude estimate; to compute it accurately at 
level (p/R) 4 one needs to include the two-instanton contribution to the effective operator. 
Note that the results for and agree well with estimates from QCD sum rules 
24], BR E7TI; our value for d^ is consistent with estimates based on measurements of the 



structure function g 2 [f|. 

Summary. The predictions of the instanton vacuum for nucleon matrix elements rele- 
vant to structure functions can be summarized as follows: 

• At twist-2 level, the quark and antiquark distributions are of order unity in the 
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packing fraction. In this order they can be computed in the effective chiral theory 
without including instanton contributions to the twist— 2 operators, i.e., replacing 
covariant by short derivatives in eq. ([!]), or dropping the path-ordered exponential 
in the light-cone operator, eq.(Q). The twist-2 gluon distribution is of order (p/i?) 4 ; 
to compute it one needs to take into account at least two-instanton contributions to 
the effective operators. 

• The instanton vacuum implies a hierarchy of twists: Large — that is, (p/R)° - 
contributions are found in operators of lowest spin (twist-4, in our case), while the 
contributions to operators of higher spin (twist-3 and 2) are suppressed. The reason 
for this pattern is the 0(4)-symmetry of the single instanton. 

The instanton vacuum provides a consistent framework for describing the non-perturbative 
input necessary for a complete understanding of DIS experiments. The key element, 
which makes possible a systematic approach to non-perturbative phenomena, is the small 
parameter p/R inherent in this picture. The "quarks-antiquarks only" approximation 
for twist-2 operators, in connection with the chiral quark soliton model of the nucleon, 
gives a very successful description of the twist-2 quark and antiquark distributions of 
the nucleon, both polarized and unpolarized 0. As to higher twists, one may hope that 
increasing accuracy of the measurements of polarized and unpolarized structure functions 
(power corrections) or, possibly, semi-inclusive measurements, will allow to test the specific 
predictions of the instanton vacuum more accurately. 

Much remains to be done on the theoretical side. In particular, one should refine 
the effective operator approach to be able to compute also parametrically small matrix 
elements, first of all the twist-2 gluon distribution. In addition to taking into account 
the two-instanton contributions to the effective operator this requires to compute also 
the effective quark action to higher orders in Mp, since many properties — for example, 
the correct realization of the QCD equations of motion — depend on the consistency of 
the definitions of effective operators and the effective action. Work in this direction is in 
progress. 
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